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Abstract 



We study the off-equilibrium response and correlation functions and the corre- 
sponding fluctuation-dissipation ratio for a purely dissipative relaxation of an 
O(N) symmetric vector model (Model A) below its upper critical dimension. 
The scaling behavior of these quantities is analyzed and the associated univer- 
sal functions are determined at first order in e = 4 — d in the high-temperature 
phase and at criticality. A non trivial limit of the fluctuation-dissipation ratio 

lf - tN + 2 ) + o<S). 



is found in the aging regime X°° = — I 1 
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I. INTRODUCTION 



The time evolution of a system relaxing towards equilibrium is characterized by two 
different regimes: a transient behavior with off-equilibrium evolution, for t < £ eq , and a 
stationary equilibrium evolution for t > t cq . In the former a dependence of the behavior 
of the system on initial conditions is expected, while in the latter homogeneity of time and 
time reversal symmetry (at least in absence of external fields) are recovered; dynamics of 
fluctuations are thus described in terms of "equilibrium" dynamics, with a characteristic 
time scale diverging at the critical point (critical slowing down) 0. 

Consider a ferromagnetic model in a disordered state for the initial time t = 0, and quench 
it at its critical temperature. During the relaxation a small external field h is applied at 
x = after a waiting time s. At time t the order parameter response to h is given by the 
response function R^(t,s) = 5(cj) x (t))/5h(s), where is the order parameter and (■) stands 
for the mean over stochastic dynamics. Since the system does not reach the equilibrium this 
function will depend both on s (the "age" of the system) and t. This behavior is usually 
referred to as aging and was firstly noted in spin glass systems [0. 

To characterize the distance from equilibrium of an aging system, evolving at a fixed 
temperature T, the fluctuation-dissipation ratio (FDR) is usually introduced f|,|4]]: 

TR^t,s) 

X ^ S) d s C x (t,s)' {lA) 

where C x (t, s) = (0 x (t)0o(s)), the two-time correlation function. When the waiting time 
s is greater than t eq the dynamics is homogeneous in time and the fluctuation-dissipation 
theorem leads to X x (t, s) = 1. This is no longer true in the aging regime ||. 



In recent years, several works P-[T3H have been devoted to the study of the FDR for 
systems exhibiting domain growth flI3| , or for aging systems such as glasses and spin glasses, 
showing that in the low-temperature phase X(t, s) turns out to be a non-trivial function of 
its two arguments. In particular analytical and numerical studies indicate that the limit 

X- = lim hmX x=0 (t, S ), (1.2) 

vanishes throughout the low-temperature phase both for spin glass and simple ferromagnetic 
system J^ffl. 

Only recently P,[l5|-|I9[| attention has been paid to the FDR, for non-equilibrium, non- 



disordered, and unfrustrated systems at criticality. It has been argued that the FDR ( |1.2|) 
is a novel universal quantity of non-equilibrium critical dynamics. Correlation and response 
functions were exactly computed in the simple cases of a random walk, a free Gaussian field, 
and a two-dimensional XY model at zero temperature and the value X£L = 1/2 was found 
@. The same problem has been addressed for the d-dimensional spherical model |17|], for 



the one dimensional Ising-Glauber chain JT^JTH| and Monte Carlo simulation was done for 



the two dimensional Ising model ||17|| . In all cases X^L has values ranging between and 



while for some urn models a different range has been found |2U||. Also the scaling form for 



R x=0 (t,s) was rigorously established using conformal invariance |f2~T| . 

In this work we investigate the non-equilibrium correlation and response functions and 
the associated FDR for the O(N) ferromagnetic model with purely dissipative relaxation 
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dynamics (Model A of Ref. JT|) both at the critical point and in the high-temperature 
phase, using a field-theoretical approach (never applied so far — see concluding remarks of 



Ref. ||17||), at first order in an e-expansion. 

The paper is organized as follows. In Section [TI| we briefly introduce the model, the 
scaling forms and the Gaussian (mean-field) result. In Section |IT| we derive the e-expansion 
for the FDR for all values of s and t. Finally in Section [IV] we summarize our results and 
discuss some points needed further investigation. In Appendix we report some useful details 
on the one loop calculation. 



II. THE MODEL 



Let us consider the purely dissipative relaxation dynamics of a iV-component field y?(x, t) 
described by the stochastic Langevin equation (Model A of Ref. [|TJ) 

^(x,*) = -o/^ + £(x,0, (2.1) 
where T~C[<p] is the Landau-Ginzburg Hamiltonian 

WM = / d d x ^-(d V f + \W + , (2.2) 

Q the kinetic coefficient, and £(x, t) a zero-mean stochastic Gaussian noise with 

(&(x,*)f J -(x',0) = 2^5(x - x')5(t - t'^ij. (2.3) 



The equilibrium correlation functions, generated by the Langevin equation ( [2.1| ) and 
averaged over the noise £, can be obtained by means of the field-theoretical action [^,|5|] 



S[<p,0\ 



dt / d d x 



-dip sn[ip] 

ip— + Utp—i ipilip 

at d(p 



(2.4) 



where <£>(x, t) is the response field. 

In Ref. |23| this formalism was extended in order to incorporate a macroscopic initial 
condition into Eq. (|2.4| ): one has also to average over the initial configuration v 2 o( x ) = <p{t = 
0, x) with a weight e~ Ho ^ given by 



H [<po] = |d d x|(^ (x)-a(x)) 2 . 



(2.5) 



This specifies an initial state a(x) with correlations proportional to r . In this way all 



response and correlation functions may be obtained, following standard methods p2| , p5| , by 
a perturbative expansion of the functional weight e - ^^'^ 4 "^ ^ ^. 

The propagators (Gaussian two point correlation and response functions) of the resulting 
theory are 



(<&(q, 8)<pj(-q, t)) = 5 i:j B?(t, s) = 6 tj 9(t-s)G(t-s) 



(<pi(q, s)v?j(-q, t)) = 6ijC°(t, s) 



q 2 + r 



G(\t-s\)+fe^--lj G(t 



(2.6) 
(2.7) 



where 



G(t) 



-Q(<j 2 +r )t 



(2.8) 



It has also been shown that r 1 is irrelevant (in the renormalization group sense) for large 



times behavior 23 



A. Scaling forms 



When a ferromagnetic system is quenched from a disordered initial state to its critical 
point, the correlation length grows as t 1 ^, where z is the dynamical critical exponents [|TJ. 
So in momentum space, applying standard scaling arguments, all the universal functions 
depend only on the two products q z t and q z s. 

In particular we expect the scaling forms 



R q (t, s) = q- 2+ ^ +z F R (Slq z (t - s),t/s) 



CJt,s) 



-2+?? 



F c (nq z (t-s),t/s), 



(2.9) 
(2.10) 



where 9 is the initial-slip exponent of response function related to the initial-slip exponent 



of the magnetization 8' and to the autocorrelation exponent A c |24| by the relation p3 

9' = 9 + z~ 1 (2 - z - rf) = z-\d - A c ). (2.11) 

The functions Fr(i/,x) and Fc(y,x) are universal apart the normalizations for small 
arguments. These functions are regular functions of both arguments, and for large x they 
behave as: 



F R (y,x) 
Fc(y,x) 



F™(y) + 0(x- 1 ) 

F?(y) 



x 



+ 0(x~ 



(2.12) 
(2.13) 



so that, for s — >• 0, these scaling forms reduce to ones of Ref. p3[ . We would also mention 
that transforming Eqs. fl2.9|J2.10|) in the real x space (to this end we have to assume that 
Fr and Fc are rapidly decreasing functions of y for y — > oo) and setting x = one easily 



obtain the scaling forms reported in Ref. ||T7| , pT| . They also reduce to the equilibrium ones 
II when t ~ s > 1. 



Let us introduce 26 



JlRq(t,s) 
<9 s C q (t, S) 



(2.14) 



then, from the scaling forms above is simple to show that, assuming F R (y,x) = 0(y), Wx, 
the FDR may be written as 



lim lim X Q= n{t. s) — (1 

-'— >oo t— >oo M ' 



-i gr(o) 



(2.15) 
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B. Gaussian FDR 

For the Gaussian model we know exactly the response and correlations functions, so we 
can evaluate the FDR (in |§ the related quantity X x has been considered, see Sec. |TV| ). 
From Eqs. ( |2.6| , |2.7p and definition ( |2.14j ) we have 

X;(t,s) = 1 = (l + e-^^ + ^r^e- 2 ^ 2 ^^ 1 . (2.16) 

If the theory is off-critical (r ^ 0) the limit of this ratio for s — > oo is 1 for all values 
of q, in agreement with the idea that in the high-temperature phase all modes have a finite 
equilibration time, so that equilibrium is recovered and as a consequence the fluctuation-dis- 
sipation theorem holds. For the critical theory, i.e. ro oc T — T c = 0, if q ^ the limit ratio is 
again equal to one, whereas for q = we have X® =Q (t, s) = 1/2. This analysis clearly shows 
that the only mode characterized by aging, i.e. that "does not relax" to the equilibrium, is 
the zero mode in the critical limit. 



III. ONE-LOOP FDR 

The aim of this section is the computation of the non-equilibrium response and correlation 
functions for the purely dissipative dynamics of the iV-vector model at one-loop order. We 
use here the method of renormalized field theory in the minimal subtraction scheme. The 
breaking of homogeneity in time gives rise to some technical problems in the renormalization 



procedure in terms of one-particle irreducible correlation functions (see |23| and references 
therein) so our computation is done in terms of connected functions. 



At one-loop order we have to evaluate, taking also into account causality p2 |, the three 
Feynman diagrams in Figure p], one for the response function and two for the correlation 
one. In terms of these diagrams we have 

R q (t, s) = R% s) - ^r^go(a) + O(g 2 ) , 

C q (t, s) = C% s) - ^ 9o ((b) + (c)) + 0{gl) . (3.1) 

In order to evaluate the FDR at criticality we have to set in this perturbative expansion 
ro = (massless theory). We also set Tq 1 = 0, since it is an irrelevant variable [23], and 



Q — 1 to lighten the notations. The first step in the calculation of the diagrams is the 
evaluations of the critical "bubble" B c (t), i.e. their common one-particle irreducible part. 
We have, in generic dimension d, 

f d d q „ , , 1 {2t) 1 - d ' 2 F(d/2-l) , rf/ , , , 

BJt) = / — -4-r CJt, t) = -— \ ; W/Q = -N d v — U 1 -^ 2 , 3.2 

v ; J (27r) d qK ' ; d/2-1 (4n) d / 2 2 d l 2 ' v ; 

2 

where N d = - — -77- — — — . Note that the equilibrium contribution to BJt) is zero for d > 2. 
(4:ii) d ' 2 T(d/2) 

Let us consider t > s in the following. We may write 
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(c) 

FIG. 1. Feynman diagrams contributing to the one-loop response (a) and correlation func- 
tion ((b) + (c)). Response functions are drawn as wavy-normal lines, whereas correlators are normal 
lines. A wavy line is attached to the response field and a normal one to the order parameter. 



(a) = j o dt'R° q (t,t')B(t')R° q (t',s), 

(b) = J Q dt'R q (t,t')B(t')C° q (t',s), 

(c) = j o dt'R° q (s,t')B(t')C° q (t',t), 



(3.3) 



where we set r = in R q and C q . 

Performing the integration and expanding in powers of e we find for the response function 



R q (t, s) = G{t -s)(l + 9o^^ In I) + 0(e\ g 2 ), 



(3.4) 



and for the correlation function 

™ / G(t - s) - G(t + s) / „ N + 2 , t 
C 9 (t, s) = 1 (l + go— In - 



_ N + 2G(t + s) 2 2 _ 2 

K -^f(2q 2 s) + O(e 2 ,g 2 ), 



9o 
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where 



/(«) 



/"df ln£e ? + (1 -e")lnt; 
Jo 



(3.5) 



(3.6) 
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and go = N^go- Note that /(0) = 0, /'(0) = —2 and /(f) has the following asymptotic 
expansion, for ?j>1 



e v ( 1 2 fc! ' 

/ v = -2 - 1 + - + - + ... + - + .. . 

V \ V V V 



(3.7) 



In order to obtain the critical functions we have to set the renormalized coupling equal 
to its fixed point value. At first order in e EH 



9o = 9 



6 



iV + 8 



0(e 2 ) 



(3.8) 



Finally we get (called P 



TV 



TV+2 > 
TV+8' 



R q {t, s) = G{t - s) ( 1 + In t ) + 0(e 2 



(3.9) 



G(t-s) -G(t + s) 
9 2 



Pi 



TV 



1 + e— In - 

4 s 



n p N G(t + s) 



4 g 2 



/(2g 2 S )+0(e 2 ), (3.10) 



that are fully compatible with the scaling form given in the previous section, with 

F R (y,x) = e-y + 0(e 2 ), (3.11) 

and 

2y 



Fc(y,x) = e 



i + 4/ 



x — 1 



e -y^ + 0(e 2 ) 



(3.12) 



In particular we recognize the exponent 9 = + 0(e 2 ) in agreement with Ref. P3" 

z = 2 + 0(e 2 ), i] = 0(e 2 ) as expected, and that Fn(y, x) is not affected by 0(e) corrections. 
It is also easy to find that 



TV 



+ 0(e 2 ) . 



(3.13) 



Computing the derivative with respect to s of the two-time correlation function and 
taking its ratio with the response function we have 



X-\s) = l+e 



2r 



Pve 



-2? 



q 2 s 



f(2q 2 s) + 2f(2q 2 s) 



+ 0(e 2 ). 



(3.14) 



Note that, at least at this order, the result is independent of the observation time t. Using 
the large v behavior of f(v), cf. Eq. (|3.7| ), we find that the limit of the FDR for s — >• oo is 
equal to 1 for all q ^ 0. Instead for q = we have (using (|3.6|) ) 



P 



^o = -Jl-^)+0(e 2 ) 



(3.15) 



in agreement with Eq. ( |2.15| ) and with the scaling forms ( P-llj ) and (|3.13|) . 
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Taking into account the effect of the mass r (deviation from critical temperature) in 
the previous computations, one obtains for the non-critical bubble (contributing to the mass 
renormalization) 



Bit) = N d 



^_-ir( d /2)r(l- d /2,2r ot ) 



dj2-\ 



(3.16) 



where T(x,y) is the incomplete T function | 27[] . Using this expression it is possible to 
determine, as previously done, correlation and response functions. We report the basic 
formulas in the Appendix. The final result is obtained computing the ratio X q in terms of 
the renormalized parameters of the theory. It is then trivial, but algebraically cumbersome, 
to show that is equal to 1 for all q in the high temperature phase. 



IV. DISCUSSION 

In this work we considered the off-equilibrium properties of the purely dissipative relax- 
ational dynamics of an iV-vector model in the framework of field theoretical e-expansion. 



We computed at first order in e the FDR, as defined in ( 2.14 ) as a function of the waiting 
time s and of the observation time t both at criticality and in the high-temperature phase. 
The main result is that the ratio is always 1 unless at criticality for q = 0, when it takes 
the value 

fl ~ = if 1 _£^ +0( ^ ) . (4.1) 
q "° 2 V 4 iV + 8 ' 



To compare our result with some particular limit considered in the literature |3|,|17 



we 



have to relate this quantity to the analog in the real x space. The following heuristic 
argument may be useful to realize that the two ratios are exactly equal, i.e. 

XZo = X?=o ■ (4-2) 

We may rewrite the FDR in real x space as a mean value of that in momentum space with 
a weight given by _R q : 

= J-dV.C q (M) = Jd'gi^M)^ 
x=0 " Tfd*qR*(t,s) ~ JdiqR^s) " ^ " 1 " ] 

Now, since we expect i? q oc e~ q ^~ s \ in the limit s,t — > oo (in the right order) X~l will 
take contributions only for the q = mode, i.e. apart a normalization, the weight function 
-R q is a S(q). However we note that at the first order in e the equality (|4.2j) is identically 
satisfied by our result since i2q oc e - q2{t - s) , cf. Eq. 

In the limit N — > oo Eq. (|47T| ) reduces to X°° = 1/2 — e/8 + 0(e 2 ) that is the same as 
the expansion of the result for the spherical model near four dimension |17[| . 

A Monte Carlo simulation of the two dimensional Ising model gave for the FDR the value 
X°° = 0.26(1) [|l]J, qualitatively in agreement with our result for e = 2, X°° = 5/12 < 1/2. 
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To have a reliable quantitative prediction the knowledge of higher loop contributions is 
required. 

Setting e = 1 for N — 1, one obtains 11/24. This number, that is a rough estimate of 
the actual three-dimensional value, may be measured in numerical or experimental works. 

This work may be easily extended to more realistic models than those previously consid- 
ered in literature, contributing to the understanding of out-of-equilibrium dynamic phenom- 
ena, currently under intensive investigation, by means of the powerful tools of perturbative 
field theory. 
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APPENDIX: DETAILS OF COMPUTATIONS 

We summarize here the main analytical results useful for the computation of correlation 
and response functions for r > at one loop. Again one has to perform all the needed 
integrations over the times, as in Eq. ( |3.3| ) with the free field correlator and response func- 
tion ( |2.6| , p.7|) . At variance with the critical theory a renormalization of the parameter r 
is now required to cancel dimensional poles both in R q and C q . 

Let us introduce the function 

Y(t)= fdrB{r) , (Al) 
Jo 



and 



W{t) = C &t G{-2t) B(t) , (A2) 
Jo 



where G(t) and B(t) are given in Eq. ( |2.8| ) and ( 3.16 ), respectively. In terms of Y and W 
we obtain (for t > s and Tq -1 = 0, O = 1) 

(a) = G(t-s)[Y(t)-Y(s)\, (A3) 

and 

(b) + (c) = ^— {G{t - s)[Y(t) - Y(s)} - G(t + s)[Y(t) + Y(s)} + 2G(t + s)W(s)} . 
1 + r o 



(A4) 



In the following with Y and W we mean also their analytic continuation in d. 
An explicit computation leads to 
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dl 2-2 

r ' 



Y(t) = {(2r t + d/2 - l)[r(l - d/2) - r(l - d/2, 2r Q t)] + {2r Q t)^e-^} 



(AS) 



and 
W(t) 



d/2-l 



, )( |jr)f/ , (j2 , ro {G(-2t)[r(l - d/2) - r(l - d/2, 2r t)] - (g 2 /^/ 2 - 1 A(l - d/2, 2g 2 t)} , 

(A6) 



/ 



where we introduced 



(for v < its analytic continuation has to be considered). 
Expanding ( |A5 , [A6 ) in e = 4 — d, we obtain 



(A7) 



-2r t 



2(4n f 2 Y{t) = —(2r t + 1) - (2r t + l)[ 7 (2r t) - lnr ] + 1 + + °( e ) 

6 aTqi 



(A8) 



and 



2(47r) d /V + r )^(t) 



= _-[r G(-2t) + g 2 ] + g 2 [lng 2 - 5{2q 2 t)} - r Q G{-2t)[j(2r t) - lnr ] + 0(e) 



where 



and 



j(v) = 1 + e~ v (lnv + -) + jTdf lnf e"* 



It is easy to find that f(v) in Eq. (|3.6|) is related to 8{v) by 



/(«) = 2 



1 + In w — 5(v) — 



(A9) 



(MO) 



(All) 



(A12) 



Plugging Eqs. (|A8|) and ( |A9| ) into Eq. ( |A3| ) and (|A4]) and then into Eqs. (|3.1] ) it is easy to 
realize that to cancel dimensional poles both in R q (t,s) and C q (t,s) a renormalization of 
the bare mass r is sufficient (at least in the case Tq 1 = we are considering) 



r = Z r r with Z r = l + N ^ 2 ^ d/2 ] - + 0{gl) 



(A13) 



in agreement with what one would expect from the corresponding static field-theory (see, for 
instance, Ref. |p5|). All the previously stated results easily follow from explicit expressions 
given above. 
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